function [ Fc ] = tyre_fc()
a1=-22.1;
a2=1011;
a3=1078;
a4=1.82;
a5=0.208;
a6=0;
a7=-0.354;
a8=0.707;
a9=0.028;
a10=0;
a11=14.8;
a12=0.022;

syms alfa Fz
Y=0;
C=1.3;
D=a1*Fz^2+a2*Fz;
B=((a3*sin(pi*a4*atan(a5*Fz)/180))/(C*D))*(1-a12*abs(Y));
E=a6*Fz^2+a7*Fz+a8;
dSh=a9*Y;
dSv=(a10*Fz^2+a11*Fz)*Y;
phi=(1-E)*(alfa*180/pi+dSh)+(E/B)*180*atan(B*(alfa*180/pi+dSh))/pi;

Fc=simplify(D*sin(C*atan(B*phi))+dSv);
end

